AD-A008  280 


AN  INTERACTIVE  SYSTEM  FOR  LONG-PERIOD 
SEISMIC  PROCESSING.  TECHNICAL  REPORT 
NO.  11.  VELA  NETWORK  EVALUATION  AND 
AUTOMATIC  PROCESSING  RESEARCH 

Frode  Ringdal,  et  al 

Texas  Instruments,  Incorporated 


Prepared  for: 

Air  Force  Technical  Applications  Center 
Advanced  Research  Projects  Agency 


December  1974 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


UNCLASSIFI 


J'CUH.7  V CL  r iC  A VlQ»  Or  Tw;u  >L  rwhmn  Del# 

I RFPORT  DOCUMENTATION  PAGE 


[t  n:im»FS 


I*  TlTLE  (enO  *•*(•  le) 


|1  GOVT  ACCCIMOM  I 


READ  INSTR» TTlONS 
BEfOWr.  ■~nWPl.ETl>  c FORM 


» t»»i  or  mfcrt  t frmoti  ccvtsiD 


AN  INTERACTIVE  SYSTEM  FOR  LONG-PERIOP 
SEISMIC  PROCESSING  U 


TecKnical 


|7 


Frode  Ringdal  and  Jeffrey  S.  Shaub 


t PERFORMING  ORG  REPORT  NUMBER 

ALEX(Ol)-  TR-74-11 

t CONTRACT  OR  GRANT  MUMICRi 

F08606-74-C-0033 


* I'tn-ORMINr.  OR'iASlI  ATIOM  NAME  AMD  AODREJS 

Texas  Instruments  Incorporated 
Equipment  Group 

Dallas.  Texas  75222 

< i Controlling  office  name  ano  address 

Advanced  Research  Projects  Agency 

Nuclear  Monitoring  Research  Office 


10.  PROGRAM  ELEMENT  PROJECT  TASK 
AREA  A WORK  uN  T NUPRERS 


VELA  T/4705/B/ETR 

12  REPORT  DATE 

10  December  1974_ 

is.  HUMBER  Or  PAGES  S 


4.  I'CM  ^*7|H.*s  I A DD*  fSVif  trom  Control  tog  OWf)  I IS.  SECURITY  CLASS  (ol  (Me 


Air  Force  Technical  Applications  Center 
VELA  Seismolgical  Center 
Alexandria,  Virginia  22314 


UNCLASSIFIED 


TSp.  DCCL  ASSIFiCATION'  DOWNGRADING 
SCHEDULE 


|is  DISTRIBUT'D*  ST  • T«-«ENT  (el  ihle  Report) 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED 


IT.  DISTRIBUTION  ST  ATEMEMT  (el  Hie  etetrecl  entered  In  Block  /0.  II  Plllotont  Bom  Report) 

Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

• US  D,p,rtm,nl  ol  Common, 

Spnngfiold,  VA.  22151 

'•  SlIRRI.CMCPTARr  NCTFS 

ARPA  Order  No.  2551 


•5  KEY  WORCf  fC>r*Hru»  on  revere*  ei«f#  if  rereeeery  •»</  Idenlffy  by  Week  w^kerj 

Seismology  Interactive  processing 

Signal  analysis  Computer  graphics 

Matched  filtering 
Very-long  period  experiment 

20  / GSTRkCT  fCenl'nus  ->n  re-eree  eld#  U n#ce«#ery  mn<S  frfenflfy  by  btock  nuefctf)  XhiS  TCpOTt  tlCSCTibcS  2.11 

interactive  graphics  processing  system  developed  for  the  PDP-15  computer 
at  the  Seismic  Data  Analysis  Center.  The  system  is  designed  to  process 
long-period  seismic  data  and  includes  options  for  interactive  bandpass 
filtering,  matched  filtering,  spectral  analysis  and  parametei  measurements. 
The  main  intention  of  the  program  development  has  'ween  to  investigate  the 
| feasibility  of  utilizing  interactive  graphics  for  event  detection  purposes  in 

DD  , r;;:MT,  1«73  e0,7,0N  1 HOV  ” '*  obeolete  , UNCLASSIFIED  

I SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dele  Knlere.h 

PHCB  SUBJECT  TO  CHANGE 


- t ■ 


1^1 


a 


ItCUAlTY  CLAWiPlCATIO*  O'  TMtl  P 


20.  Continued 

a potential  world-wide  seismic  surveillance  system.  An  evaluation 
cf  .he  expected  operational  capabilities  of  the  system  has  been  carried 
out,  using  data  recorded  by  the  Very  Long-Period  Experiment  network 
as  a source.  Conclusions  from  this  evaluation  and  recommendations 
for  future  study  are  presented. 


UNCLASSIFIED 


••CURlTV  CL  AU4PIC  ATlON  Or  THI|  PAOCPtltwi  Data  mlar+f) 


f 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED 

ALEX(01)-TR-74-11 


D 

D 

0 

il 

0 

D 

D 

0 

0 

0 

D 

0 

0 

|] 

0 

D 


AN  INTERACTIVE  SYSTEM  FOR  LONG-PERIOD  SEISMIC  PROCESSING 


TECHN.CAL  REPORT  NO.  11 


VELA  NETWORK  EVALUATION  AND  AUTOMATIC  PROCESSING  RESEARCH 

Prepared  by 

F rode  Ringdal  and  Jeffrey  S.  Shaub 


TEXAS  INSTRUMENTS  INCORPORATED 
Equipment  Group 
Post  Office  Box  6015 
Dallas,  Texas  75222 


Prepared  for 

AIR  FORCE  TECHNICAL  APPLICATIONS  CENTER 
Alexandria,  Virginia  22314 


Sponsored  by 

ADVANCED  RESEARCH  PROJECTS  AGENCY 
Nuclear  Monitoring  Research  Office 
ARPA  Program  Code  No.  4F10 
ARPA  Order  No.  2551 


10  December  1974 


Acknowledgment:  This  research  was  supported  by  the  Advanced 

Research  Projects  Agency,  Nuclear  Monitoring  Research  Office, 
under  Project  VELA -UNIFORM , and  accomplished  under  the  tech- 
nical direction  of  the  Air  Force  Technical  Applications  Center 
under  Contract  No.  F08606 -74 -C -00 33. 


[ I 

. 


1 


Equipment  Group 


I ID 

Id 


ABSTRACT 


This  report  describes  an  interactive  graphics  process- 
ing system  developed  for  the  PDP-15  computer  at  the  Seismic  Data  Analysis 
Center.  The  system  is  designed  to  process  long-period  seismic  data  and 
includes  options  for  interactive  bandpass  filtering,  matched  filtering, 
spectral  analysis  and  parameter  measurements.  The  main  intention  of  the 
program  development  has  been  to  investigate  the  feasibility  of  utilizing 
interactive  graphics  for  event  detection  purposes  in  a potential  world-wide 
seismic  surveillance  system.  An  evaluation  of  the  expected  operational 
capabilities  of  the  system  has  been  carried  out,  using  data  recorded  by  the 
Very  Long-Period  Experiment  network  as  a source.  Conclusions  from  this 


evaluation  and  recommendations  for  future  study  ere  presented. 


Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  re  ^onsible  for  inform?  tion  contained 
herein  which  has  been  supplied  by  other  organizations  or  ccnt'actors,  and  this 
document  is  subject  to  later  revision  as  may  be  necessary.  Tha  views  and  con- 
clusions presented  are  those  of  the  author?'  and  should  not  be  interpreted  as 
necessarily  representing  the  official  policies,  either  expressed  or  iixplied,  of 
the  Advanced  Research  Projects  Agency,  the  Air  Force  Technical  Applications 
Center,  or  the  US  Government. 
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SECTION  I 
INTRODUCTION 


This  report  describes  an  interactive  processing  system 
developed  by  Texas  Instruments  Incorporated  for  the  purpose  of  analyzing 
long-period  seismic  signals.  The  system  utilizes  the  interactive  graphics 
facilities  of  the  PDP-15/50  computer  at  the  Seismic  Data  Analysis  Center 
(SDAC)  in  Alexandria,  Virginia.  The  program  development  was  conducted 
as  part  of  the  System  Study  task  of  the  VELA  Network  Evaluation  and  Auto- 
matic Processing  Recearch  program. 

The  main  intention  of  the  program  system  has  been  to  investi- 
gate the  feasibility  of  utilizing  interactive  graphics  for  event  detection 
purposes  in  a potential  world-wide  seismic  surveillance  system.  Although 
the  software  package  is  general  in  nature,  an  interface  to  one  specific 
seismic  system  (the  Very  Long-Period  Experiment  (VLPF)  network)  has  been 
designed,  and  a preliminary  evaluation  has  been  conducted  using  data 
recorded  by  this  network  as  a source. 

Interactive  processing  is  today  a very  important  aspect  of 
numerous  computer  applications.  It  provides  an  efficient  means  for  a user 
to  comprehend  his  data  base,  to  direct  a computer  in  its  operations  upon 
that  data  base,  and  to  examine  the  results  of  those  operations  - all  within 
an  appropriate  time  interval.  The  principal  advantages  of  interactive 
processing  are: 

• It  reduces  the  waiting  time  between  intermediate  processing 

steps,  thus  increasing  productivity 

• It  reduci  s the  need  for  hard-copy  output  because  a video 

di  lay  of  intermediate  results  is  sufficient  in  many  applications 
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• It  provides  an  efficient  means  to  retain  human  judgment 

in  the  analysis  loop,  and  thus  avoids  the  problems  inherent 
in  fully  automating  analytical  decisions. 

Interactive  processing  is  particularly  well  suited  for  those 
applications  that  are  characterized  as  a series  of  sub-processes  with  active 
intermediate  decision  points.  Seismic  signal  analysis  belongs  to  this  class 
of  problems.  Typical  intermediate  decision  points  are  exemplified  as 
follows: 

• Data  quality  control;  elimination  or  correction  of  bad  data 
segments. 

• Alignment  of  signal  traces  for  beamforming. 

• Selection  of  the  "best"  bandpass  filter  or  matched  filter 
from  a filter  library. 

• Selection  of  a signal  peak  for  magnitude  measurements. 

• Selection  of  time  windows  for  computing  quantities  such 
as  seismic  noise  level  and  the  AR  and  AL  discriminants 
(Brune  et  al.  , 1963). 

• Rapid  visual  control  of  detection/no-detection  decisions  on 
individual  signal  traces. 

In  addition,  several  non-routine  seismic  signal  processing 
techniques  may  benefit  greatly  from  interactive  processing.  Examples  in- 
clude the  complex  Cepstrum  technique,  identification  of  later  phases  (such 
as  pP),  and  detection  association  techniques  for  network  processing.  For 
a discussion  of  these  and  related  topics  we  refer  to  Sax  (1974) 
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The  interactive  system  describes  in  this  report  deals 
primarily  with  standard  processing  techniques  such  as  bandpass  filtering, 
linear  chirp  or  reference  waveform  matched  filtering,  computatior  of 
power  spectra  and  measurements  of  selected  event  parameters.  In  Section 
II  we  present  junctional  description  of  the  system.  Section  III  gives 
some  examples  of  application  and  includes  a brief  evaluation  of  the  pro- 
cessing effectiveness  in  a potential  surveillance  mode.  Finally,  conclu- 
sions an.'  recommendations  for  future  study  are  presented  in  Section  IV. 

A documentation  of  the  developed  software  has  been  issued 
separately  from  this  report  (Ringdal  and  Shaub,  1974).  This  documentation 
also  contains  a user  description  of  the  programs,  including  a step-by-step 
solution  of  a sample  problem. 


1-3 


SECTION  II 

DESCRIPTION  OF  THE  ILPPS  INTERACTIVE  SYSTEM 


A.  PROGRAM  PURPOSE 


The  overall  purpose  of  the  Interactive  Long-Period  Processing 
System  (ILPPS)  is  to  provide  an  interactive  graphics  capability  for  detecting 
and  analyzing  seismic  surface  waves.  Although  the  software  has  been  design- 
ed to  operate  on  a specific  computer  system  (the  PDP-15/50),  the  design  phi- 
losophy is  sufficiently  general  to  apply  to  other  potential  configurations.  The 
ILPPS  system  is  primarily  intended  to  support  two  objectives: 

• To  provide  an  analyst  with  a convenient  tool  to  perform  inter- 
active fignal  analysis  for  research  purposes. 

• To  establish  an  operationally  flexible  system  which  might  be 
used  to  process  a large  number  of  events  in  a world-wide  seis- 
mic surveillance  network. 


In  order  to  achieve  these  objectives,  the  ILPPS  system  has  been 
designed  with  emphasis  on  establishing  a convenient  sequential  processing  ca- 
pability to  accomodate  the  requirements  of  an  operational  surveillance  system. 
At  the  same  time,  enough  flexibility  has  been  built  in  to  allow  the  analyst  to 
select  non-standard  processing  functions  at  various  intermediate  steps  in  the 
analysis. 

B.  SYSTEM  CONFIGURATION 
1.  Hardware 


i 


The  ILPPS  system  has  been  implemented  for  the  PDP-15/50 
computer  system  located  at  SDAC.  lhis  configuration  comprises  the  folio 
ing  hardware  features: 


w- 


Central  Processing  Unit  with  64  K words  memor  r (24  K are 
used  by  ILPPS). 


Standard  peripheral  equipment  including  four  magnetic  tape 
stations,  card  reader  and  line  printer. 


Two  disk  drives;  one  large  capacity  disk  capable  cf  storing 
10  million  words  and  one  fixed-head  disk  of  256  K word 
capacity. 


A video  display  interactive  graphics  console  with  keyboard, 
light  pen,  and  interrupt  pushbuttons. 

A VARIAN  electrostatic  printer/plotter  and  a CALCOMP 
digital  plotter. 


For  a more  detailed  description  of  the  system  hardware  con- 
figuration we  refer  to  Teledyne/Geotech  (1973). 


2. 


Software 


The  initial  implementation  of  ILPPS  runs  under  the  RSX 
Plus  II  operating  system  for  the  PDP-15  computer.  RSX  provides  a multi- 
programming environment  and  thus  allows  for  the  execution  of  other  tasks 
concurrently  with  ILPPS.  This  capability  is  clearly  of  considerable  impor- 
tance when  operating  in  an  interactive  mode,  since  the  computer  will  be 
inactive  between  processing  steps,  while  the  analyst  makes  his  decisions. 
This  wait  time  will  for  most  applications  be  considerably  larger  than  the 
actual  CPU  time  utilized. 
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The  ILPPS  software  is  coded  primarily  in  FORTRAN.  The 
only  exceptions  are  certain  input/output  functions  for  which  assembly 
language  has  been  utilized. 

3.  Data  Base 

In  the  initial  version  of  the  ILPPS  package,  the  capability 
exists  to  process  dat?  recorded  by  the  Very  Long-Period  Experiment  (VLPE) 
Network.  Table  II- 1 lists  the  individual  VLPE  stations.  For  further  descrip- 
tion of  the  network,  we  refer  to  Lambert  et  al.  (1973). 

C.  FUNCTIONAL  DESCRIPTION 

1.  General  Features 

The  basic  structure  of  the  ILPPS  software  package  is  out- 
lined in  Figure  II - 1 . The  main  features  of  the  system  may  be  described  as 
follows: 

• An  interactive  graphics  console  which  serves  as  the  physical 
interface  between  the  computer  system  and  the  analyst. 

• A system  supervisor  which  provides  a command  language  to 
control  processing  module  execution  within  the  interactive 
environment. 

• A set  of  independent  processing  modules  which  perform  analysis 

functions  such  as  the  selection,  display,  and  filtering  of  seismic 

waveforms. 

This  modular  design  approach  provides  for  a simple  logical 
structure  of  the  system.  It  furthermore  facilitiates  future  extensions  of  the 
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TAELE  II- 1 

VERY  LONG-PERIOD  EXPERIMENT  (VLPE) 
STATION  AND  LOCATION 


Station  Name  Designator 


Charters  Towers,  Australia 
Chiang  Mai,  Thailand 
Fairbanks,  Alaska 
Toledo,  Spain  TLO 

Eilat,  Israel  EIL 

Kongsberg,  Norway  KON 

Ogdensburg,  New  Jersey  OGD 

Kipapa,  Hawaii  KIP 

Albuquerque,  New  Mexico  ALQ 

La  P"  z,  Rolivia  ZLP 

Matsushiro,  Japan  MAT 


Latitude 


20.09S 
18.  79N 
64.  90N 
39.  86N 
29.  55N 
59.65N 
41.07N 
21.42N 
34.  94N 
16. 50S 
36.  54N 


Longitude 


146. 

26E 

98. 

98E 

148. 

01W 

4. 

02W 

34. 

95E 

74. 

62W 

158. 

02W 

106. 

46  W 

68. 

13W 

138. 

21E 

Off-Line 

Program 


Fixed 

Head 

Disk 


PSCAN 


FILTER 


Display 
Waveforms 
From  Disk 


Perform 

Signal 

Analysis 


Select 
Event  And 
Time  Window 


'Interac- 
tive  Graph 
ics  Displav 


SUPERVISOR 

System 

Linkage 


TI  INTERACTIVE  LONG-PERIOD  PROCESSING  SYSTEM 

GFNFPAL  FLOW 


[ 
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system  functions,  either  by  riodifying  existing  processes  or  by  including  addi- 
tional software  modules.  This  is  especially  imoortant  in  view  of  the  experi- 
mental nature  of  the  ILPPS  system. 

A number  of  additional  design  considerations  apply  to  the 
ILPPS  system.  The  most  important  of  these  are  listed  below: 

• The  data  base  is  completely  disk-oriented;  thus  event  data 
are  stored  cn  disk  by  an  off-line  program  prior  to  system 
execution.  The  capability  exists  to  input  new  event  data  to 
the  disk  while  existing  events  are  being  processed. 

• Although  the  system  is  initially  designed  to  process  VLPE 
network  data,  the  off-line  approach  to  system  input  makes 
adaption  to  o.her  types  of  seismic  data  easy. 

• Data  commumcation  between  distinct  processing  modules  is 
performed  via  temporary  disk  files,  in  order  to  minimize 
memory  requirements. 

• For  ease  of  operation,  analyst  control  of  the  program  flow 
is  maintained  mainly  via  pushbutton  interrupts.  At  each 
decision  point,  a function  selection  panel  is  displayed;  giving 
key-words  to  explain  to  the  analyst  the  available  options. 

• All  keyboard  entries  may  be  made  in  free  format.  A compre- 
hensive error  checking  system  recognizes  possible  syntax 
errors  as  well  as  input  values  that  are  out  of  range. 

It  might  be  no.ed  that  the  ILPPS  system  in  most  practical 
applications  will  accept  as  input  long-period  seismic  waveforms  corres- 
ponding to  a known  event.  Normally  the  event  will  have  been  detected  by 
a short-period  station  or  network;  thus  an  approximate  location  and  origin 
time  are  known.  This  approach  is  reasonable,  since  short-period  P-wave 
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detection  is  usually  much  easiei  than  identification  of  surface  waves.  How- 
ever, the  ILPPS  programs  are  general  enough  to  process  waveforms  for 
which  no  event  information  is  provided,  in  case  this  should  be  desired. 
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2.  Processing  Options 

The  following  paragraphs  detail  some  of  the  processing 
options  within  the  ILPPS  system.  A further  description  of  some  examples 
of  application  will  be  presented  in  Section  III. 

a.  System  Supervisor 

• Creates,  updates,  and  retains  user  procedures.  Previously 
created  procedures  may  be  accessed  from  a disk  storage 
file. 

• Accesses  system  processing  modules  either  individually 
or  in  a procedural  mode. 

• Maintains  a complets  record  of  system  functions  performed, 
with  hard-copy  output  to  the  printer. 

The  internal  structure  and  design  of  the  ILPPS  system 
supervisor  was  influenced  substantially  by  that  of  the  Numerical  Analysis 
Problem  Solving  System  (Roman,  197  3). 

b.  DPSCAN  Processing  Module 

• Lists  summary  of  all  events  in  current  data  base 

• Provides  random  access  to  any  event  in  the  data  base 

• Allows  rapid  scanring  of  individual  events  by  displaying  one 
station  component  at  a time  (maximum  4096  seconds  time 
window). 
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c.  SELLV  Processing  Module 


Provide.,  random  access  to  any  individual  data  channel  (i.  e.  , 


the  vertical,  transverse,  or  radial  component  trace  for  a 


specified  station  and  event) 


Allows  analyst  option  to  select  a time  window  for  processing 


by  use  of  a moving  cursor 


Stores  selected  traces  on  a disk  file  for  later  access  by  the 


IILTER  routine 


Allows  simultaneous  selection  of  any  number  and  combination 


of  traces,  stations  and  events;  the  only  restriction  being  the 
amount  of  computer  core  storage  available  for  the  display  file. 


d.  FILTER  Processing  Module 


Displays  simultaneously  all  traces  selected  by  the  SELEV  routine 


Performs  bandpass  filtering,  linear  chirp  matched  filtering, 
or  reference  waveform  matched  filtering  of  any  individual 


trace 


Includes  a sequential  processing  capability  for  selecting 
optimum  chirp  filter  ’ength 


Can  adjust  scaling  par  meters  (horizontal  or  vertical)  and 


select  various  levels  of  annotation 


Provides  a save-restore  capability  to  and  from  disk  for 


the  displayed  picture 


Displays  signal  or  noise  power  spectra  and  chirp  impulse 


response. 
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• Includes  an  option  to  apply  a set  of  narrow-band  filters  to  a 

wave-form  to  determine  its  dispersion  characteristics. 

• Measures,  interactively,  the  following  parameters: 

- M (Love  or  Rayleigh  waves) 

s 

- RMS  Noise  levrel 

- AR  or  AL  parameters  (Brune  et  al.  , 1963) 

• Retains  measur  d parameters  in  a disk  storage  file. 

3.  Program  Execution 

All  individual  analysis  programs  are  executed  under  control 
of  the  system  supervisor.  The  analyst  may  initiate  any  one  of  these  by 
entering  a PERFORM  statement  while  the  system  is  in  supervisory  state. 

The  system  supervisor  also  includes  options  to  create  user  procedures 
combining  various  sequences  of  calls  to  system  modules  and  to  execute 
such  procedures  by  i single  request.  Figure  II - 2 shows  an  example  of  the 
system  log  in  a simple  case  where  the  analyst  has  requested  the  DPSCAN, 
SELEV  and  FILTER  modules  individ  ually. 

Within  each  analysis  package,  program  control  is  maintained 
via  the  pushbutton  interrupt  panel.  This  provides  lor  rapid  selection  and 
execution  of  processing  optiins  and  minimizes  the  necessary  manual  inter- 
action. As  an  aid  to  the  analyst,  a set  of  key  words  describing  each  of 
six  possible  interrupt  options  is  displayed  at  each  decision  point. 

Figure  II- 3 shows,  as  an  example,  the  key  word  system  within 
the  FILTER  program  module.  This  figure  illustrates  the  hierarchy  used 
in  ILPPS  to  organize  the  various  program  options.  For  example,  in  order 
to  filter  the  current  work  trace  with  a narrow  band  filter  (0.040-0.050  Hz), 
the  analyst  hits  pushbuttons  2,  1 and  4 in  that  sequence  (from  panels  2.  , 

2.  1 and  2.  1.  1 in  the  figure,  respectively). 


TI-SAPSS  LOGGED  ON  t <| ✓ 1 9/ AT  13:58:29  , ENTER  COMMAND, 


PERFORM  CPSCAN 

SYSTEM  FUNCTION  OPSCAN  INITIATED 
EVENT  LIST  REOUESTED,  ( Y , N ) ? 

ENTER  DESIRED  EVENT  NUMBER 

....  13 

....  .PERFORM  SELEV 
SYSTEM  FUNCTION  SELEV  INITIATED 
EN1ER  EVENT  SEQUENCE  NUMBER  FROM  LIST 

....  13 

ENTER  STATION  NUMBER  (1,15) 

ENTER  COMPONENT  NUMBER  (1,3) 

ENTER  NUMBER  OF  COPIES  TO  BE  SAVED  j TS  = 1337  TL  = 1322  SEC 

ENTER  EVEfT  SEQUENCE  NUMBER  FROM  LIST 

*.!!’  .PERFORM  FILTER 
SYSTEM  FUNCTION  FILTER  INITIATED 
TRACE  3 AMP=  146U.5  NM  PERs  23,2  SEC  MS=  a. 93 

PERFORM  LOGOFF 


TI-SAPSS  LOGGED  OFF  » 4/19/74  AT  11*  2:30 


FIGURE  II -2 

EXAMPLE  OF  SYSTEM  LOG  OUTPUT 
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The  net  vffect  of  this  approach  to  the  function  selection  problem 
is  to  have  the  pushbutton  panel  act  as  a set  of  programmable  function  keys. 

The  hierarchical  organization  of  the  processing  options  makes  it  feasible 
to  access  easily  a large  number  of  functions  within  this  framework. 
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SECTION  III 

OPERATIONAL  ASPECTS 

A.  EXAMPLES  OF  APPLICATION 

In  this  section,  we  present  a number  of  examples  to  illustrate 
some  of  the  capabilities  of  the  ILPPS  system.  These  examples  are  not  intend- 
ed to  describe  how  to  operate  the  system  - for  this  we  refer  to  the  software 
documentation  (Ringdal  and  Shaub,  1974).  However,  they  should  give  some  in- 
sight into  data  selection,  simultaneous  display  capabilities  and  signal  analysis 
options. 

Figure  III-- 1 shows  the  first  page  of  the  event  directory  listing 
the  event  data  sets  that  are  currently  available  on  the  disk.  Epicentrai  infor- 
mation has  been  obtained  from  short-pt  riod  data  (e.  g.  , the  LASA  or  NORSAR 
seismic  bulletins).  Station  codes  indicate  (by  a 0 in  the  appropriate  column) 
which  of  15  VLPE  stations  have  data  available  for  the  particular  event. 
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An  example  of  an  event  waveform  record  for  station  nu  Tiber  9 
(ALQ)  (vertical  component,  unfiltered)  for  an  event,  from  central  Asia  is  dis- 
played in  Figure  III-2.  The  event  annotation  (two  top  lines)  has  been  derived 
from  short-period  information,  and  estimated  arrival  times  of  long-period  P, 

S,  LQ,  and  LR  waves  are  marked.  The  tick  marks  on  the  time  axis  occur  at 
100-second  intervals;  the  total  time  period  covered  is  4096  seconds.  The 
traces  are  scaled  automatically,  with  the  scale  factor  displayed  at  the  bottom 
of  the  screen.  The  right-hand  part  of  the  screen  shows  the  six  presently  avail- 
able processing  options  corresponding  to  the  six  pushbuttons  of  the  display  unit. 
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FIGURE  III-  1 

DISPLAY  OF  ILPPS  EVENT  DIRECTORY 
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DISPLAY  OF  A SAMPLE  WAVEFORM 
(4096  SECONDS) 


III  - 3 


In  Figure  III- 3,  a simultaneous  display  of  four  different 
station  recordings  of  the  same  event  is  shown.  Each  station  is  represented 
by  its  vertical  component.  The  time  window  in  each  case  is  1024  seconds. 

The  annotation  below  the  traces  identifies  event  name,  station  number, 
component  number,  station  azimuth  and  distance  Irom  epicenter,  data 
start  time,  filter  band,  and  scale. 

Figure  III-4  gives  a : imultaneous  display  of  the  three  rotated 
components  (LRV,  LQT,  LRR,  respectively)  of  station  number  11  (MAT)  for 
a central  Asian  ev  nt.  In  this  case,  all  tra:es  have  been  bandpass  filtered; 
the  vertical  trace  by  a 0.  020-0.060  Hz  filter;  the  two  horizontal  traces  by 
a 0.015-0.060  Hz  filter. 

F.gure  III- 5 shows  an  application  of  a linear  chirp  filter  to 
a central  Asian  event  recorded  by  station  6 (KON,  vertical  component).  The 
top  trace  is  the  unfiltered  waveform;  the  middle  trace  is  the  same  waveform 
filtered  with  a narrow-band  filter  (0.040-0.050  Hz)  and  the  third  trace  is 
obtained  by  applying  a linear  chirp  of  600  seconds  over  a frequency  band  of 
0.020-0.060  Hz  to  the  original  trace. 

Figure  IiI-6  is  identical  to  Figure  III-5,  except  that  the 
impulse  response  of  the  chirp  filter  (time-inverted)  is  displayed  in  the  middle 
trace. 

Figure  III-7  shows  a waveform  recorded  by  station  9 (ALQ) 
ir.  the  top  trace,  and  the  log  power  spectrum  of  the  same  event  in  the  bottom 
trace.  Frequency  increments  are  linear,  with  the  rightmost  point  corres- 
ponding to  the  Nyquist  frequency  of  0.  25  Hz. 

Figure  IH-8  illustrates  reference  waveform  filtering.  The  top  tra 
is  the  vertical  component  of  station  8 (KIP)  for  a presumed  explosion  from 
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FIGURE  III -4 

DISPLAY  OF  THREE  ROTATED  COMPONENTS 
(FROM  TOP  TO  BOTTOM:  LR  V,  LQT,  LRR) 
OF  A VLPE  STATION 
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FIGURE  III  - 5 


EXAMPLE  OF  AN  ORIGINAL  EVENT  WAVEFORM 
(TOP  TRACE),  FILTERED  BY  A NARROW-BAND  FILTER 
(MIDDLE  TRACE)  AND  BY  A LINEAR 
CHIRP  FILTER  (BOTTOM  TRACE) 
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F IGURE  III  - v 

DISPLAY  OF’  A SEISMIC  WAVEFORM 
AND  ITS  POWER  SPECTRUM 
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FIGURE  III-  8 

DISPLAY  OF  A RECORDED  WAVEFORM  (TOP  TRACE) 

A REFERENCE  EVENT  FROM  THE  SAME  LOCATION  (BOTTOM  TRACE) 
AND  THE  CROSS-CORRELATION  BETWEEN  THE  TWO  (MIDDLE  TRACE) 
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eastern  Kazakh;  the  bottom  trace  is  the  reference  waveform;  which  is  a 
much  larger  event  from  the  same  site.  The  middle  trace  is  the  cross 
correlation  between  the  top  and  bottom  traces.  It  is  observed  that  the 
signal-to-noise  ratio  of  the  middle  trace  is  substantially  improved  over 
that  of  the  top  trace. 

Figure  III  - 9 pictures  the  outputs  of  a suite  of  narrow-band 
filters  applied  to  a selected  recording.  The  multiple  filter  option  is  designed 
to  give  an  efficient  way  of  computing  group  velocity  curves  for  various  trans- 
mission paths.  Figure  III- 10  shows  the  second  part  of  this  option,  the 
Hilbert  transforms  of  the  same  traces  with  cursor  marks  (set  by  the 
analyst)  at  the  local  maxima.  Finally,  in  Figure  III -11  the  travel  times 
are  translated  into  group  velocities.  The  box  surrounding  a point  indicates 
how  the  analyst  may  select  points  on  the  actual  group  velocity  curve  for 
the  path  in  question. 

B.  EVALUATION 

1.  Description  of  Experiment 

The  ILPPS  system  might  conveivably  be  used  in  the  future 
either  for  research  purposes  or  in  a potential  surveillance  mode.  Clearly, 
both  applications  demand  a ligh  quality  of  the  data  analysis.  Apart  from 
this,  the  requirements  differ.  In  a research  application,  the  most  valuable 
feature  is  flexibility  and  adaptability  to  non-standard  operations.  In  a 
surveillance  system,  the  real-time  aspect  implies  that  processing  efficien- 
cy is  of  primary  importance. 

We  describe  here  an  experiment  that  was  conducted  in  order 
to  evaluate  the  effectiveness  of  ILPPS  in  a potential  world-wide  surveillance 
system.  Data  recorded  by  the  VLPE  network  was  utilized  for  the  evaluation; 
the  stations  of  this  network  were  previously  listed  in  Table  II - 1 . 
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FIGURE  III  - 10 

A SET  OF  ENVELOPES  (HILBERT  TRANSFORMS) 

OF  A SUITE  OF  NARROW-BAND  FILTFRED  TRACES. 

THE  CURSOR  MARKS  INDICATE  ANALYST  TIME  DICKS 
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FIGURE  III- 11 

CONSTRUCTION'  OF  A GROUP  VELOCITY 
CURVE  VIA  HIE  INTERACTIVE  DISPLAY 
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A total  of  10  events  ^specified  in  Table  III  - 1 ) were  processed. 
All  of  these  events  occurred  during  January  of  1973.  Only  6 VLPE  stations 
were  operating  reliably  during  this  period  (stations  number  2,  5,  6,  8,  9 
and  11);  consequently,  the  processing  was  restricted  to  those  stations.  Only 
the  vertical  component  of  each  station  was  analyzed  in  this  study. 


In  order  to  conduct  a realistic  experiment,  we  included  in 
the  processing  all  the  options  of  ILPPS  that  could  reasonably  be  expected 
to  be  used  in  an  actual  surveillance  system.  Typically,  each  station  wave- 


form was  processed  as  follows: 


Extract  the  station  waveform  from  the  Event  Data  File  on  Disk. 


Select  a time  windov  for  processing  that  includes  the  pre- 
dicted arrival  of  the  Rayleigh  wave  train  in  a 2.  5-3.  9 kr.i  sec 
velocity  window.  For  stations  of  epicentral  distance  less  than 
60  degrees,  a time  window  of  1024  seconds  is  usually  selected, 


otherwise,  a 2048  second  time  window  is  used. 


Perform  bandpass  filtering  with  a 0.  020-0.  060  Hz  filter, 


Perform  linear  chirp  filtering  with  chirp  window  lengths 
selected  by  the  analyst  using  the  bandpass  above.  Typically, 
5-10  different  chirp  lengths  are  tried. 


Perform  reference  waveform  filtering  using  a fixed  and 


previously  selected  event  as  a reference.  (See  Table  III  - 1) 


Compute  magnitude  on  the  bandpass  filtered  trace.  (If 
no  signal  was  detected,  this  will  be  a "noise  magnitude"  which 
is  useful  as  an  upper  bound  on  the  actual  station  magnitude. ) 
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TABLE  III- 1 

LIST  OF  PROCESSED  EVENTS 


***  processed  Events'  *** 


event  *t>  ’ 

■i£ 

OP  I G DATE 

TIKE 

LAT 

LON 

ME 

STATION  CODES 

LX+CEKAP+ 

11 

1 / ^ 3/ 7 3 

15.05.1 8 

39N 

72E 

9.8 

101000100001111 

L X + C E n A r.  + 

12 

1/03/73 

17.09.93 

39N 

72E 

3.8 

101000*00001111 

LX+CE’  A\t 

37 

1/12/73 

17.35.99 

38N 

70Z 

3.7 

101000100001111 

IX+CESAp* 

55 

1/18/73 

06.97.28 

39N 

99E 

9.6 

1010001000011 1 1 

LX+CENAP+ 

59 

1/19/73 

15.10.02 

32N 

6BE 

5.0 

101000100001111 

Lx  + CE*»a\+ 

61' 

1/19/73 

16.92.91 

35N 

7 1 E 

3.6 

101000100001111 

L >.  ♦ CT  •'  A + 

65 

1/20/73 

19.31.59 

3 1 N 

67E 

9,0 

101000100001111 

LX  + Ctf.  ALt 

o7 

1/21/73 

03.23.53 

91  N 

7 1 E 

9.3 

101000100001111 

Lx  + CE  AF  + 

71 

1/23/73 

11.31.98 

90N 

9 1 E 

9.Q 

101000100001111 

LX  + CE  • AP  + 

75 

1/29/73 

03.20.20 

9 1 N 

82E 

5.1 

1010001000011 1 1 

***  REFERENCE 

EVENTS 

* * ★ 

EVENT  '.t  ■ E ORIG  DATE 

TIME 

LAT 

LON 

MB 

STATION  CODES 

LXtCEl  AP+  6 1/02/73 

IX  + Cf  ai>+  K 1/03/73 

22.25.57 

19.31.09 

3 1 N 
39N 

88E 
7 IE- 

5.2 

5,5_ 

1010001000011 1 1 
101000100001111 

Explanation: 
Station  codes: 
(Stations  1 to  15): 
Orig  Date/Time: 
LA  1 -LON: 

EV#: 
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0 : Data  available 

1 : No  data  available 

Event  origin  date  and  time  ; GMT 
Epicentral  latitude  and  longitude  (degrees) 

Event  number  of  processed  events  (1-10) 

Ml  is  a master  event  used  for  Station  6 

M2  is  a master  event  used  for  Stations  .2,  5,  8,  9»  11 
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• If  a signal  was  detected,  compute  the  AR  discriminant  on 
the  bandpass  filtered  trace  (Brune  et  al.  , 1963). 

• Compute  RMS  noise  and  signal-to-noise  ratio  for: 

- tine  unfiltered  trace 

- the  bandpass  filtered  trace 

- the  "best"  chirp  filtered  trace  (only  if  the  event  was  detected) 

- the  master  filtered  trace  (only  if  the  event  was  detected). 


• If  a detection/no -detection  decision  is  questionable,  go  back 
to  process  the  full  event  trace  (4096  seconds)  to  get  a better 
visual  indication  of  noise  versus  signal  characteristics. 

• Store  event  parameters  on  disk  together  with  data  quality 
and  detection  indicators  set  by  the  analyst. 

r 

2.  Results 

Table  III  - 2 gives  a breakdown  of  the  processing  times  in 
minutes  for  the  station  data  analyzed.  The  time  required  to  process  one 
component  has  a relatively  wide  range  (from  3 to  8 minutes  if  faulty  data 
channels  are  disregarded).  The  most  typical  values  are  3 or  4 minutes  for 
one  trace,  with  an  average  processing  time  of  4 minutes  for  good  channels. 

As  anticipated,  the  longest  time  was  spent  in  processing  components  with 
marginal  detections. 

In  Table  III  - 3 we  present  a breakdown  of  the  time  required 
to  process  the  individual  subtasks  that  together  form  the  complete  pro- 
cessing of  one  vveform.  These  times  are  all  estimated  "typical"  times, 
and  add  up  to  about  3 to  7 minutes  for  best  and  worst  case,  respectively. 

The  following  observations  may  be  made: 
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TABLE  III- 3 

BREAKDOWN  BY  PROCESSING  SUBTASK  OF  TIME 
REQUIRED  TO  ANALYZE  A TYPICAL  WAVEFORM 


Type  of  Operation 

Total  Time 
(including 
CPU  time) 

Estimated 

CPU 

Time 

Select;on  of  waveform  segment 

20-30  sec 

5 sec 

Bandpass  filtering  (512  pts  or  1024  pts) 

10-20  sec 

3-6  sec 

Chirp  filtering  (scanning  of  5-10  chirp 
lengths) 

1-3  min 

15-45  sec 

Reference  waveform  filtering 

20-30  sec 

10-15  sec 

Parameter  measurements 

30-60  sec 

5-10  sec 

Selection  and  filtering  of  "long  trace" 
(not  usually  performed) 

60  sec 

20  sec 

'Total  Time"  includes: 


CPU  time 


Computer  wait  time  (e.  g.  , for  disk  data  retrieval) 


Time  required  for  analyst  entry  of  numbers  or  commands 


A reasonable  (although  short)  time  for  the  analyst  to  judge 
the  results. 
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• CPU  time  (i.  e.  , time  used  for  actual  computations  by  the 
Central  Processing  Unit  of  the  PDP-15  computer)  constitutes 
approximately  25  percent  of  the  total  time  actually  spent  on 
processing  an  event. 

• The  single  subtask  requiring  the  most  time  is  the  chirp 
filter  scanning  (1-3  minutes). 

• Both  the  selection  of  waveform  segments  (20-30  seconds) 
and  the  parameter  measurements  (30-60  seconds)  consume  a 
substantial  amount  of  time  especially  relative  to  the  corres- 
ponding CPU  time  utilized,  which  is  on  the  order  of  5 seconds 
in  each  case. 

It  is  clearly  of  considerable  interest  to  examine  possible 
ways  to  speed  up  routine  processing.  This  problem  will  be  addressed  in 
detail  in  Subsection  III-C. 

We  proceed  to  present  the  actual  data  analysis  results  from 
the  evaluation  experiment.  Since  our  primary  purpose  was  not  to  evaluate 
the  VLPE  network  as  such  or  the  filtering  techniques  applied;  but  rather 
the  interactive  analysis  concept,  the  actual  processing  results  are  of 
secondary  interest.  However,  we  still  include  these  results  for  completeness 
and  to  verify  that  the  interactive  procedure  produces  results  of  acceptable 
quality. 

Table  III -4  presents  the  effects  on  detection/no -detection 
decisions  when  examining  matched  filter  output  traces  instead  of  conventional 
bandpass-filtered  data.  As  expected,  the  matched  filter  makes  only  a mar- 
ginal difference  in  the  number  of  detections.  Out  of  the  34  station-events 
that  were  not  detected  with  bandpass  filtering,  two  had  a chirp  filter  detection. 
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One  additional  event  improved  from  a marginal  detection  to  a clear  detection 
when  applying  a chirp  filter.  The  results  of  master  waveform  filtering 
were  slightly  inferior  to  the  chirp  results;  probably  due  to  the  application 
of  a fixed  master  waveform  in  all  cases  rather  than  trying  to  determine  an 
optimum  reference  trace. 

We  would  like  to  point  out  that  our  criteria  for  declaring 
detection  were  very  conservative  in  this  evaluation;  i.  e.  , we  required  a 
clearly  visible  dispersion  pattern  in  order  to  accept  a detection  on  band- 
passed  data,  and  a matched  filter  detection  was  only  declared  if  a clear 
peak  was  observed  in  the  expected  signal  arrival  window. 

Matched  filter  detection  performance  in  this  study  is  some- 
what inferior  to  that  determined  by  Strauss  (1974)  using  a more  extensive 
data  base  of  central  Asian  events  recorded  by  the  VLPE  network.  We 
attribute  this  difference  to  the  slightly  less  strict  detection  criteria  applied 
in  his  study. 

The  gains  in  signal-to-noise  ratio  (SNR)  for  the  matched 
filtered  waveforms  relative  to  the  bandpassed  traces  are  shown  in  Table 
HI _ 5 for  all  waveforms  with  a detection.  Note  that  the  SNR  values  of  the 
bandpassed  traces  represent  maximum  zero -to -peak  value  of  the  signal 
relative  to  the  root-mean-square  (RMS)  value  of  the  noise;  thus  in  a sense 
they  are  biased  high.  However,  the  chirp  and  master  filtering  gains  are 
true  values,  based  upon  compatible  processing  in  each  case. 

The  following  observations  may  be  made: 

• Chirp  filter  gains  vary  from  less  than  1 dB  to  more  than 

6 dB,  with  an  average  value  of  2.  7 dB. 
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TABLE  III -5 

MATCHED  FILTERING  GAINS 
FOR  DETECTED  EVENTS 


Event 

Number 

Station 

Number 

Distance 

(degrees) 

SNR 

(bandpass) 

dB 

Chirp 

Gain 

dB 

Master 

Gain 

dB 

Chirp 

Length 

Seconds 

I 

2 

30.9 

8.  8 

2.  2 

2.  9 

280 

4 

2 

16.  1 

19.  5 

1.6 

3.  8 

75 

5 

2 

30.  8 

16.  1 

1.  5 

2.  7 

280 

9 

2 

22.  7 

15.  1 

0.  5 

1.4 

100 

10 

2 

26.4 

29.0 

5.0 

2.  0 

230 

6 

5 

30.  8 

14.  1 

0.  6 

0.  5 

200 

I 

6 

43.  5 

9-  6 

2.  8 

3.0 

120 

9 

6 

51.9 

14.7 

0.  8 

1.  6 

320 

10 

6 

47.  ? 

23.  1 

2.  1 

6.  8 

450 

4 

8 

87.4 

9.  3 

4.  5 

0.  7 

900 

5 

8 

110.  1 

19.  8 

5.  9 

2.  2 

790 

10 

8 

96.  3 

17.  4 

3.  5 

4.  5 

900 

4 

9 

106.0 

11.  4 

2.  9 

3.  3 

1630 

5 

9 

112.  2 

23.  2 

6.4 

5.  5 

1600 

I 

II 

51.  5 

9.  6 

2.  1 

- 

200 

9 

11 

36.7 

10.  4 

2.  7 

3.0 

250 

10 

11 

43.  1 

18.  9 

0.  7 

-1.6 

4 50 

2.  7 dB 

3.0  dB  (excluding  station  11) 
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• Master  waveform  filtering  gains  have  about  the  same  range 
as  the  chirp  gains,  except  for  station  11,  which  had  a master 
waveform  with  a poor  signal-to-noise  ratio.  The  average 
gain  for  master  waveform  filtering,  not  including  station  11, 
is  3.0  dB. 

• As  expected,  matched  filter  gains  are  generally  highest  for 
the  most  distant  stations  (which,  of  course,  usually  have  the 
most  dispersed  waveforms.) 

• In  several  cases,  there  was  a significant  dif f c rence  (up  to 

+ 5 dB)  between  master  filtering  gain  and  chirp  filtering  gain. 

Table  III  - 5 also  shows  that  the  length  of  the  best  linear  filter 
varies  considerably  even  within  narrow  epicentral  distance  ranges.  This 
implies  that  it  would  be  difficult  to  predict  the  optimum  chirp  filter  length 
for  a given  event,  therefore,  fairly  extensive  calibration  data  would  be 
necessary  in  order  to  cut  down  on  the  number  of  iterations  necessary  to 
determine  the  best  linear  chirp. 

The  Rayleigh  wave  magnitudes  of  the  processed  events  are 
listed  in  Table  III -6 . In  those  cases  where  no  detection  was  declared,  the 
number  listed  indicates  the  "noise  magnitude,  " which  is  an  upper  bound 
for  the  actual  event  magnitude  at  the  station.  The  "network  magnitudes,  " 
averaged  over  all  detecting  stations  for  each  event,  are  also  included  in 
the  table. 

The  low  noise  magnitudes  of  station  5 (EIL)  and  the  poor 
detection  performance  of  this  station  indicate  a discrepancy  n the  amplitude 
response  of  the  vertical  component  of  this  station  for  the  processed  events. 
This  assertion  is  supported  by  the  fact  that  the  radial  component  of  this 
station  was  found  to  have  more  detections  and  substantially  higher  mag- 
nitudes than  the  vertical  component. 
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C.  DISCUSSION 

The  evaluation  experiment  described  in  the  preceding  subsection 
demonstrates  that  the  ILPPS  system  provides  processing  results  of  a quality 
comparable  to  those  obtained  by  conventional  methods.  In  addition,  interactive 
processing  of  an  event  may  be  completed  in  a considerably  shorter  time  span 
than  is  possible  with  batch  processing,  and  with  significantly  improved  conven- 
ience for  the  analyst. 


However,  the  time  factor  still  remains  the  most  fundamental 
question  when  one  considers  the  application  of  ILPPS  in  a potential  seismic 
surveillance  system.  At  an  average  processing  time  of  4 minutes  per  station, 
about  11/2  hours  would  be  required  per  event  in  a hypothetical  25  station  net- 
work, and  considerably  more  time  would  be  expended  if  processing  of  both 
Love  and  Rayleigh  waves  were  required.  These  requirements  seem  prohibi- 
tive for  routine  analysis  of  all  detected  events,  but  might  be  acceptable  for 
processing  events  cf  special  interest  in  a surveillance  system. 

In  this  context,  it  is  important  to  remember  that  the  processing 
times  presented  in  the  ILPPS  evaluation  are  based  upon  event  analysis  in  a 
fully  interactive  mode.  By  this  we  mean  that  the  analyst  has  specifically  se- 
lected each  option  and  each  parameter  setting  in  each  particular  processing 
case.  Clearly,  it  would  be  possible  to  obtain  a significant  reduction  in  time 
requirements  by  adapting  a more  automated  analysis  procedure  with  analysis 
intervention  only  at  a few  specific  decision  points.  An  example  of  such  a semi- 
automated  processing  system  could  be  as  follows: 

• Initially,  the  full  waveform  (2  048  points)  is  displayed,  and  the 

analyst  selects  either  a "manual"  processing  mode  (as  described 
previously)  or  an  "automated"  mode,  which  is  described  in  the 
following. 
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The  expected  signal  arrival  window  is  extracted,  and  the  fol- 
lowing traces  are  displayed  automatically. 

The  original  trace 

The  trace  filtered  with  a standard  bandpass  filter 
The  trace  filtered  with  a default  chirp  filter  (based  upon 
regional  information) 


The  trace  correlated  with  a reference  waveform  taken 


from  a library. 


The  analyst  then  has  the  option  to  request  additional  processing 
but  presumably  in  most  cases  he  will  make  a detection/no  de- 
tection decision  based  upon  the  displayed  data. 


A parameter  measurement  routine  is  then  invoked.  All  mea- 


surements are  automated,  but  with  an  analyst  override  capability. 


After  each  component  is  processed,  control  passes  automatically 
to  the  next  data  channel  by  direct  disk  access. 


A procedure  as  outlined  above  would  probably  make  it  feasible 
to  reduce  processing  time  to  less  than  1 minute  in  the  'standard'  cases  (i.  e.  , 
if  no  analyst  override  action  occurs).  Although  this  type  of  standard  process- 
ing might  be  somewhat  inferior  to  the  fully  interactive  procedure  (since  no  it- 
eration on  chirp  filter  lengths  would  normally  be  performed),  it  might  still  be 
adequate  in  most  cases.  Clearly,  the  full  interactive  processing  machinery 
would  always  be  available  to  ensure  optimum  processing  of  events  of  special 


interest. 


The  semi-automated  capability  outlined  above  is  not  available 
in  the  initial  implementation  of  the  ILPPS  system.  However,  the  inclusion  of 
such  an  option  would  not  require  a very  large  programming  effort,  and  should 
be  considered  in  future  applications  of  the  system. 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  puipose  of  the  ILPPS  experiment  was  to  investigate  the 
feasibility  of  using  interactive  graphics  for  processing  long-period  data 
in  an  operational  seismic  surveillance  system.  It  is  felt  that  the  following 
features  of  the  interactive  approach  have  been  definitely  demonstrated: 

• High  quality  of  results 

• Convenience  to  the  analyst 

• Minimal  intermediate  hard  copy  output 

• Short  turn-around  time  compared  to  batch  processing. 

The  one  major  question  not  fully  answered  is  whether  inter- 
active signal  analysis  is  efficient  enough  for  the  large-scale  routine  pro- 
cessing required  in  a surveillance  system.  The  average  processing  time 
for  one  station  component  during  the  ILPPS  evaluation  was  4 minutes, 

» 

including  time  for  event  selection,  bandpass  and  iterative  matched  filtering 
and  interactive  computation  of  several  event  parameters.  This  processing 
time  's  probably  prohibitive  for  routine  analysis  in  surveillance  mode. 

However,  it  is  possible  to  reduce  the  average  ILPPS  processing  time 
significantly  by  the  following  approach: 

• Establish  a semi-automated  interactive  system,  in  which 
a fairly  extensive  default  processing  may  be  performed 
automatically  if  the  analyst  so  wishes. 

• Retain  an  option  to  perform  extensive  interactive  analysis 

1 

of  difficult  cases  or  events  of  special  interest. 
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• Improve  computer  efficiency  by  various  means  (e.  g.  , 
extensive  use  of  direct  access  disk  operations  to  reduce 
wait  time). 

The  time  required  for  this  type  of  a routine  processing  could 
probably  be  reduced  to  about  1 minute  per  station  component.  This  pro- 
cessing time  would  seem  to  satisfy  the  real-time  requirements  of  a large 
scale  surveillance  system,  while  the  indicated  approach  still  retains  the 
desired  flexibility  to  perform  extensive  analysis  of  interesting  events. 

It  is  therefore  recommended  that  further  development  of 
ILPPS  be  directed  toward  establishing  a semi -automated  processing  capability 
to  supplement  the  already  existing  fully  interactive  system.  The  addition  of 
such  a capability  would  probably  not  require  a major  software  effort. 

Also,  improving  efficiency  in  the  computer  operations  and 
in  providing  more  analyst  conveniences  such  as  hard  copy  output  from  CRT 
should  be  given  high  priority. 

Finally,  severed  additional  options  may  be  included  in  the 
system  at  relatively  low  cost,  such  as: 

• The  capability  to  process  short-period  data. 

• Techniques  for  data  quality  contro.  and  spike  removal. 

• Additional  processing  techniques,  such  as  beamforming, 
three-component  processing,  complex  cepstrum  and  multi- 
channel filtering. 

• The  capability  to  interface  directly  with  a remote  seismic 
data  mass  storage,  in  the  event  such  a system  is  established. 
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The  implementation  of  some  or  all  of  the  above  options  can 
be  expected  to  provide  more  insight  into  areas  within  the  seismic  event 
detection  problem  which  are  well  suited  to  the  application  of  interactive 
graphics.  This  information  will  be  valuable  both  for  seismic  data  pro- 
cessing techniques  in  general  and  also  for  the  possible  future  operation  of 
a global  seismic  surveillance  network. 
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